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^ ■ Abstract 

1-^ I An isospin-selfconsistent approach based on the Continuum-Random-Phase- Approxi- 

' mation (CRPA) is apphed to describe the Fermi and Gamow-Teher strength distributions 

within a wide excitation-energy interval. To take into account nucleon pairing in open-shell 
nuclei, we formulate an isospin-selfconsistent version of the proton- neutron-quasiparticle- 
1/-^ I CRPA (pn-QCRPA) approach by incorporating the BCS model into the CRPA method. 

■ The isospin and configurational splittings of the Gamow-Teller giant resonance are ana- 

, lyzed in single-open-shell nuclei. The calculation results obtained for ^^^Bi, ^^Nb, and Sb 

' isotopes are compared with available experimental data. 
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1. Introduction 



O I The Fermi (F) and Gamow-Teller (GT) strength functions in medium-heavy mass nuclei 

^ I have been studied for a long time. The subject of studies of the GT-strength distribution is 
closely related to the weak probes of nuclei (single and double /5-decays, neutrino interaction, 
etc.) as well as to the direct charge-exchange reactions ((p,n), (^He,t), etc.). The weak probes 
deal with the low-energy part of the GT strength distribution (see, e.g., lU EI)- The direct 
reactions allow to study the distribution in a wide excitation-energy interval including the 
region of the GT giant resonance (GTR) (see, e.g., jSHZ] and references therein) and also the 
high-energy region (see, e.g., [Hlin])- The discovery of the isobaric analog resonances (lAR) and 
the subsequent study of their properties have allowed one to conclude on the high degree of 
the isospin conservation in medium-heavy mass nuclei. It means that the lAR exhausts almost 
the total Fermi strength and the rest is mainly exhausted by the isovector monopole giant 
resonance (IVMR) 

The Random-Phase- Approximation (RPA)-based microscopical studies of the GT strength 
distribution started from the consideration of a schematic three-level model taking the di- 
rect, core-polarization, and back-spin-flip transitions into account PT|. One or two weakly- 
collectivized GT states along with the GTR were found in the model. Realistic Continuum- 
Random-Phase-Approximation (CRPA) calculations for the GT strength distribution were per- 
formed in [12] and jT^]. Some of the states predicted in pQ] were found in [T3]. Attempts to 
describe the GT strength distribution in details within the RPA-|-Hartree-Fock model have 
been undertaken recently in ^14j. Unfortunately, the authors of fl2\ lT^ did not address the 
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following questions: 1) the GT strength distribution in the high-energy region of the isovec- 
tor spin-monopole giant resonance (IVSMR); 2) effects of the nucleon pairing; 3) the isospin 
splitting of the GTR. 

Having taken the nucleon pairing into consideration in CRPA calculations, the authors 
of [13] predicted the configurational splitting of the GTR in some nuclei and initiated the 
respective experimental search for the effect In ^3]; however, the influence of the particle- 
particle interaction in the charge-exchange channel was not taken into account. Along with the 
pairing interaction in the neutral channel, the proton-neutron interaction is taken into consid- 
eration within the proton-neutron-quasiparticle-Random-Phase-Approximation (pn-QRPA) to 
describe the double /3-decay rates (see, e.g., and references therein). Unfortunately, most 
current versions of the pn-QRPA do not treat the single-particle continuum that hinders the de- 
scription of both high-lying IVMR and IVSMR. In addition, it seems that the question whether 
the modern versions of the pn-QRPA comply with the isospin conservation has not been raised 
yet. 

The present paper is stimulated partially by the experimental results of [3HZ1 and by the 
intention to overcome (to some extent) the shortcomings of previous approaches. As a base 
we use the isospin-selfconsistent CRPA approach of [TTHT^ . where direct-decay properties of 
giant resonances have been mainly considered. We pursue the following goals: 

1) application of the isospin-selfconsistent CRPA approach to describe the F and GT strength 
distributions in closed-shell nuclei within a wide excitation-energy interval including the 
region of the isovector monopole and spin-monopole giant resonances; 

2) formulation of an approximate method to deal with the isospin splitting of the GTR; 

3) taking into account the spin-quadrupole part of the particle-hole interaction for descrip- 
tion of GT excitations; 

4) incorporation of the BCS model into the CRPA method to formulate an isospin-selfconsistent 
version of the proton-neutron-quasiparticle-Continuum-Random-Phase- Approximation (pn- 
QCRPA) approach; 

5) application of the pn-QCRPA approach to describe the F and GT strength distributions 
in single-closed-shell nuclei within a wide excitation-energy interval; 

6) examination of the GTR configurational splitting effect within the pn-QCRPA and its 
impact on the total GTR width. 

We restrict ourselves to the analysis within the particle-hole subspace, and, therefore, do not 
address the question of the influence of 2p-2h configurations (see, e.g., jTH]) as well as the 
quenching effect (see, e.g., [IHl)- However, we simulate the coupling of the GT states with 
many-quasiparticle configurations using an appropriate smearing parameter. 

The above-listed goals could be apparently achieved within the pn-QCRPA approach based 
on the density-functional method [20] • However, the authors of [201 focused their efforts mainly 
on the analysis of the low-energy part of GT strength distribution relevant for the astrophysical 
applications.^ 

The paper is organized as follows. In Section 2 the basic relationships of the isospin- 
selfconsistent CRPA approach are given. They include: description of the model Hamiltonian, 
its symmetries and the respective sum rules (subsections 2.1 and 2.2), the CRPA equations 

^We would emphasise the decisive contribution of Prof. S.A. Fayans in developing such a powerful! approach. 
We regret much on his too-early passing away. 
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with taking into account the spin-quadrupole part of the particle-hole interaction for descrip- 
tion of the GT excitations (subsection 2.3); the approximate description of the isospin splitting 
of the GTR (subsection 2.4). In Section 3 we extend the CRPA approach to incorporate the 
BCS model to describe both the nucleon pairing in neutral channels and the respective inter- 
action in the charge-exchange particle-particle channels. The generalized model Hamiltonian, 
its symmetries and respective sum rules are given in subsections 3.1 and 3.2. The standard 
pn-QCRPA equations are reformulated in terms of radial parts of the transition density and 
free two-quasiparticle propagator. That allows us to formulate the coordinate-space represen- 
tation for the inhomogeneous system of the pn-CRPA equations (subsection 3.3). As applied to 
single-open-shell nuclei, a version of the pn-QCRPA approach is formulated on the base of the 
above system (subsection 3.4). The choice of the model parameters and calculation results con- 
cerned with the F and GT strength distributions in ^°^Bi, ^°Nb, and Sb isotopes are presented 
in Section 4. Summary concerned to the approach and discussion of the calculation results are 
given in Section 5. 

2. Description of the Fermi and GT strength functions 
in closed-shell nuclei 

2.1. Model Hamiltonian 

We use a simple, and at the same time realistic, model Hamiltonian to analyze the Fermi, 
"Coulomb" (C), and GT strength functions in medium-heavy mass spherical nuclei. The Hamil- 
tonian consists of the mean field U{x) including the phenomenological isoscalar part Uo{x) along 
with the isovector Ui{x) and the Coulomb Uc{x) parts calculated consistently in the Hartree 
approximation : 

U{x) = Uo{x) + Ui{x) + Uc{x), (1) 
Uo{x) = Uo{r) + UUx); f/i(x) = ^i;(r)r(3); Uc{x) = ^Uc{r){l - r^^)). 

Here, f/o(r) and Uso{x) = Uso{r)(T\ are the central and spin-orbit parts of the isoscalar mean 
field, respectively; v{r) is the symmetry potential. The potential U{x) determines the single- 
particle levels with the energies ex {X = tt for protons and X = u for neutrons, A is the set 
of the single-particle quantum numbers, (A) = {Ij}) and the radial wave functions r~^xx{^) 
along with the radial Green's functions {rr')~^g(^x){r, r'; e). On the base of Eq. (0), one can get 
expressions relating the matrix elements of the single-particle Fermi (t'-"-') and GT (cr^r''"^) 
operators (a^ are the spherical Pauli matrices): 

^\%{u){i^^ - ^v){X^Xv) + v^v - {Uc)i^v\ = 0, (2) 

't'ff){u){i^7T - eu){X7TXu) + - {Uc)^v - ((crl)(7r) " {y)){U so) ^u} = (3) 

with iXnXi^) = IXn{r)Xu{r)dr, f^^ = J f{r)x7r{r)xAr) dr, tj^H 
v/2aTT5wH, = ^((7r)||a||(z/)). 

We choose the Landau-Migdal forces to describe the particle-hole interaction jS]. The 
explicit expression for the forces in the charge-exchange channel is: 

Hp.,, = 2j2{Fo + FiCTaCTt)TjfWi-^^6{ra - n) + h.c. (4) 

a>b 

where the intensities of the non-spin-flip and spin-flip parts of this interaction, Fq and Fi 
respectively, are the phenomenological parameters. 
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2.2. The symmetries of the model Hamiltonian and sum rules. 

The model Hamiltonian H complies with the isospin symmetry provided that 

[^,f(-)] = f/(r), (5) 

H = Ho + Hp-h, Ho = E{Ta + U{xa)), T(-)=Eri-), U^^ = EUcir,)r(-K (6) 

a a a 

Using the RPA in the coordinate representation for closed-shell nuclei, one can get according 
to Eqs. (HI), (HI)-® the well-known selfconsistency condition |221l2Sj: 

v{r) = 2Fon^~\r); n^'\r) = (0| J^rf ^^(r - r„)|0) = n''{r) - nP{r), (7) 

a 



47rr^ 



A 



where n^^\r) is the neutron excess density, are the occupation numbers {P = p,n). The 
selfconsistency condition relates the symmetry potential to the Landau-Migdal parameter Fq. 

The equation of motion for the GT operator Y^^^ = Y,ai'^fi)aT!r^ can be derived within the 
RPA analogically to Eq. © fSSJ: 

[H,Ytn = U^-\ (9) 
Ujr^ = EaU,ixa)T(-\ U,ix) = [f/.o(x),a^] + (2(Fi - Fo)n(-)(r) + Ucir))a^. (10) 

Equations (0) and ^ allow one to render some relationships being useful to check the calcu- 
lation results for the strength functions within the RPA. The strength function corresponding 
to the single-particle probing operator V^"^^ = J2a ^(^^a)'*"!^'' is defined as: 



a 
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5{?'H = ^|(s|V^W|0) 6{u;-u;,) (11) 

s 

with ujg = Eg — Eq being the excitation energy of the corresponding isobaric nucleus measured 
from the ground state of the parent nucleus. Using Eqs. (0) and pi|) one gets the relationship: 

S^^\uj)=uj-^S^(^\uj). (12) 

Here, the Fermi and "Coulomb" strength functions correspond to the probing operators T^^^ 
and tJc'\ respectively. 

The model-independent non-energy-weighted sum rule {NEWSR) for F(0"'") and GT(1+) 
excitations is well-known |21|: 

{NEWSR)j = {0\\G%^\G%^]\0) = {N - Z); Gio^=T(-); G[-^ = Y^-l (13) 



The energy-weighted sum rules (EWSR) 



{EWSR)j = (0| [H, G%>\\ |0) (14) 

are rather model-dependent and according to Eqs. (0), (jHI), and (fTUj) equal to: 

{EWSR)o = J Uc{r)n'--\r)d\ (15) 

(EWSR), = {EWSR)o - ^{0\Uso\0) + 2(Fi - Fq) J in^-\r)Yd'r (16) 



with (0|?7,„|0) = EpEx^if, - iDllill + l)(f/.o)AA. 

According to Eqs. (0), (fT^. and (fT3j). the exact (within the RPA) isospin SU(2) sym- 
metry is reahzed for the model Hamiltonian in question in the hmit Uc — Ac = {N — 
Z)~^ J Uc{r)n^~\r)d^r. In this hmit the energy and the wave function \A) of the "ideal" 
isobaric analog state (IAS) are Ea = Eq + Aq and \A) = (N — Z)^^/'^T^^^\0) , respectively. 
The "ideal" IAS exhausts 100% of {NEWSR)^. In the calculations with the use of a realistic 
Coulomb mean field, the IAS exhausts almost 100% of {NEWSR)o (with the rest exhausted 
mainly by the IVMR). Therefore, one can approximately use the isospin classification of the 
nuclear states. Redistribution of the Fermi strength is caused mainly by the Coulomb mixing 
of the IAS and the states having the "normal" isospin Tq — 1 (Tq = (A^ — Z)/2 is the isospin of 
both the parent-nucleus ground state and its analog state). In the present model the mixing 
is due to the difference Uc{r) — Ac- The Wigner SU(4) symmetry is realized in the limit 
Uso — >■ 0, Fi — i> Fo, Uc — > Ac. In this limit the energy Eg and the wave function \Gfi) of the 
"ideal" Gamow-Teller state (GTS) are Eg = Ea and = {N - Z)-^/'^Y^-^\0), respectively. 
Redistribution of the GT strength is mainly due to the spin-orbit part of the mean field. 



2.3. The strength functions within the continuum-RPA 

The distribution of the particle-hole strength can be calculated within the continuum- 
RPA making use of the full basis of the single-particle states. Being based on the above- 
described model Hamiltonian, the CRPA equations for calculations of the F (C) and GT 
strength functions can be derived using the methods of the finite Fermi-system theory |2I]. Let 
^jLS/i ~ ^a^JLSiMixa)Ta'^^ bc the multipolc charge-exchange probing operator leading to exci- 
tations with the angular moment J and parity vr = (— 1)*^. Here, VJLS^l{x) = VjL5(r)T7i5^(n), 
TjLSiJ.i''^) — ^/^YjM'm CilpsM^LM'f^sM is the irreducible spin-angular tensor operator with 
(Too = 1, o-iM = In particular, we have Vooo(^) = ^{Uc{r)), T^om = 1 and Vloi(^) = 
1, TiQi^ = cr^ for description of F(C) and GT excitations, respectively. After separation of 
the isospin and spin-angular variables, the strength functions corresponding to the probing 
operators VjIIm are determined by the following equations: 

^lU^) = Im / VjLs{r)A^f^s,JL'S'ir, r'; ^)vjl4'(^', ^) drdr', (17) 

viVs'ir.u;) = VjLs{r)6LL'6ss' + ^ E / A^7L'S',JL''S''ir,r';u)vills,,{r',uj) dr'. (18) 



Here, vjj^l{r,uj) are the effective radial probing operators, -^5=0,1 are the interaction intensities 
of Eq. (jD), {A'Krr')~^A^^j^g j^,g,{r,r';uj) are the radial free particle-hole propagators: 



,(-) . ., ,^ _ ((vr)||r,^5||(^.))((7r)||r,^,5'||(^))* 



AjLS,JL'S'irur2]uj) = 47r^ 



X 



2J+1 

{nuXu{ri)Xu{r2)9{7T){ri, + u) + n^Xn{ri)x7r{r2)9{u){ri,r2; - u)} . (19) 

The expression for A^~^^ can be obtained from Eq. (jl9p by the substitution n ^ u with 
LJ being the excitation energy of the daughter nucleus in the /J^-channel, measured from the 
ground state of the parent nucleus. For J'^ = excitations {L = S = 0) there is only one 
non-zero propagator v4gQoooo and, therefore, Eqs. (fTTj) . (fTS|) have the simplest form. Such a 
form has been used explicitly in |TB] for describing lAR and IVMR. For J'^ = 1"*" excitations 
{S = 1) propagators A^^\ are diagonal with respect to S and, therefore, Eq. (fTHj) is the 
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system of equations for effective operators Viqi and V121 . As a result, the spin-quadrupole part 
of the particle-hole interaction contributes to the formation of the GT strength function, as it 
follows from Eqs. (fT7|). (fT5|) . The use of only diagonal (with respect to L) propagators Aj^^ jn 
corresponds to the so-called "symmetric" approximation. In particular, this approximation was 
used in [T71 HH] to describe the monopole and dipole spin-flip charge-exchange excitations. As 
a rule, the use of the "symmetric" approximation leads just to small errors in calculations of 
the strength functions Sj^g in the vicinity of the respective giant resonance (see, e.g., |2SI)- 
Nevertheless, the detailed description of the low-energy part of the GT strength distribution 
appeals to the "non-symmetric" approximation. 

The F and GT strength functions Sj^\uj) (hereafter, indexes L = 0, 
S = J are omitted) calculated within the CRPA for not-too-high excitation energies (including 
the region of lAR and GTR) reveal narrow resonances corresponding to the particle-hole-type 
doorway states. Therefore, the following parametrization holds in the vicinity of each doorway 
state: 

4^^M--^Im^^%-. (20) 
vr uj - Us + 2^ s 

Here, {r^j^^)s, oOs^ and Tg are the strength, energy, and escape width of the doorway state, respec- 
tively. The values of Xj"* = Xj ^ — x^/^ = 
{IS^'\u)du - JS^j'\u)du)/{NEWSR)j, vY = y'f^ + y^^ 

J are useful to be compared with unity to check the 
quality of the calculation results. Due to the relations given in Eqs. (jllj) and (j^n|) the ratio of 
the "Coulomb" to the Fermi strengths has to be u"^ for each doorway state. 



2.4. Isospin splitting of the GT strength distribution. 

Due to the high degree of the isospin conservation in nuclei, the GT states are classified 
with the isospin Tq — 1,To,Tq + 1. In particular, Tq components of the GT strength function 
can be considered as the isobaric analog of the isovector Ml states in the respective parent 
nucleus: 

|M,To) = (2To)-^/2f(-)|Ml,To). (21) 

According to the definition (jlip . the strength function of the Tq components is proportional to 
the strength function of the isovector Ml giant resonance in the parent nucleus: 

Sir\uj) = (2To)-i5ffi(^' = ^ - Ac). (22) 

Here, S^^ji{uj') is the Ml strength function corresponding to the probing operator Mj^^ = 
Y.ai^tija'Ta^^ and depending on the excitation energy in the daughter nucleus. Within the 
CRPA the Ml strength function is calculated using the relations, which are similar to those of 
Eqs. ffT7|) . (fT5|) . and can be found, e.g., in [TTj . 

According to Eq. ()22|) . a noticeable effect of the isospin splitting of the GT strength function 
takes place only for nuclei with a not-too-large neutron excess. The suppression of the Tq and 
To -|- 1 components is the reason why the GT states \s), obtained within the RPA, are usually 
assigned the isospin To — 1, although the RPA GT states do not have a definite isospin. In fact, 
the states \s) and |M, Tq) are non-orthogonal: 

(M,To|s) = (2To)-i/2(Ml,To 

where Q^"-* is a RPA boson-type operator corresponding to the creation of a collective GT 
state |s). Therefore, one has to project |s) states onto the space of the GT states with the 



(23) 
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isospin To — 1 = T< by means of subtraction of the admixtures of To states to force the relevant 
orthogonahty condition: 

\s,n - 1) = (1 - EibMrr'^'ils) - E&m|M,To)). (24) 

M M 

This equation is vahd under the assumption that the integral relative strength a;> = Y^m^m 
of the To = T> component is small as compared to unity. In this case taking into account 
the To + 1 component is even more unimportant. We restrict the further analysis by the 
approximation that there is the only GT state in the respective RPA calculations exhausting 
100% of the NEWSR and, therefore, this state can be considered as the "ideal" GTS. Under 
these assumptions one has: 

hi, = (2To)-i/2(A/r,To|r(-)|0) =xif ; x< = 1 -x>. (25) 

Having averaged the exact nuclear Hamiltonian over the state |s,To) in the form of Eq. ()24|). 
one gets: 

00^ = X~^{XsUJs - xmujm)- (26) 
M 

According to the approximate relations and (j^Uj) . the relative strength and the energy 
of the To — 1 GTS diminish as compared with the respective RPA values. The decrease is 
determined by the zeroth and first moments of the strength function of the isovector Ml GR 
in the parent nucleus: 

x> = (AT - Z)-' f SgUo;') du'; J^^mXm = {N - Zy' f st},{uj')u' dJ . (27) 

These relations again lead to the conclusion that the value 
x> ^ 2A^^'^(N — Z)^'^ is rather small even if the value (A^ — Z) is not large. Only for nuclei 
with a minimal neutron excess (a few units) all three isospin components of the GT strength 
could have comparable strengths. 



3. The strength functions in open-shell nuclei 

3.1. The model Hamiltonian 

The interaction Hp-p in the particle-particle channels has to be included in the model Hamil- 
tonian for nuclei with open shells, in order to take the effects of nucleon pairing into considera- 
tion. We choose the interaction in the following separable form (as it is used in the BCS model) 
in both the neutral and charge-exchange particle-particle channels with the total angular mo- 
mentum and parity of the nucleon pair being J'" = 0~^, 
(L = 0,5= J): 

Hp-p = "2 ^^-^^ i^fsp'VPfBf^' ■ (28) 

Jfj. P=P,n 
P'=p,n 

Here, Gj=o,i are the intensities of the particle-particle interaction, P^^ is the annihilation 
operator for the nucleon pair: 

AA' mm' 
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where aAm(«Am) the annihilation (creation) operator of the nucleon in the state with the 
quantum numbers Am (m is the projection of the particle angular momentum). The interaction 
()28|) preserves the isospin symmetry of the model Hamiltonian. 

We use the Bogolyubov transformation to describe the nucleon pairing in the neutral chan- 
nels in terms of quasiparticle creation (annihilation) operators aAm('^Am) (see, e.g., |2H|)- As 
a result, we get the following model Hamiltonian to describe F (C) and GT excitations in the 
/3~-channel within the quasiboson version of the pn-QRPA: 

H = H, + H,.n + Hp-p, Ho= (aL)+«L -^if^P' ^^n)f^'^. (30) 

I3=p,n Am 



Here, = ^ (Ca)^ + the quasiparticle energy, — /^/9 chemical 

potential and the energy gap, respectively, which are determined from the BCS-model equations: 

iV/3 = E(2jf + l)(^Af; A, = GoE(2jf + lK^f (31) 

A A 

with v\ = |(1 — ^) and = 1 — vf. 

The total interaction Hamiltonian in both particle-hole (j3)) and particle-particle ()28|) chan- 
nels can be expressed in terms of the quasiparticle (pn)-pair creation and annihilation operators 
which obey approximately the bosonic commutation rules: 

= E (Jl^lj^rnjurn'h^m'a^m; [^2', [Ai^V] = 5^^>5^^,5jj,5^^>. (32) 

mm' 

The explicit expressions for the interactions are: 

Hp-h = E 77 E {x.X.X.'Xu'){Q'/.yQ%. (33) 

J/j ttvk'v' 

Hp-p = -J:Gj i: ix.Xu)ix.'Xu')iPiTPX', (34) 

Jfi TTUn'u' 

where Qi^^ = t[-i]^,)iu^v,Ai^, + v^uMi^y), = 45(.)(^-«-^;J^^ " '"^-^MiiV). = 

{-ly+^'Ai^'', iXnXuXn'Xw) = I XnXuXn'Xu'r"^dr. 

3.2. The symmetries of the model Hamiltonian and sum rules. 

Due to the isobaric invariance of the particle-particle interaction ()28|) . the Eq. (0) still holds 
and leads exactly to the same selfconsistency condition of Eq. ((Tj) with the only difference, that 
the proton and neutron densities are determined with account for the particle redistribution 
caused by the nucleon pairing: 

n^'\r) = ^E(2jf + iMnx'xir))'. (35) 

The direct realization of the Eq. (0) within the pn-QRPA with making use of Eqs. (j3(Jj) - (j34j) 
and T^~^ = J^iruix-n-Xu) (Q'^Z)^ leads also to the selfconsistensy condition of Eq. ((Tj) provided 
that the full basis of the single-particle states for neutron and proton subsystems is used along 
with Eq. Q for the radial overlap integrals for the proton and neutron wave functions. The 
consistent mean Coulomb field Uc{r) and {EWSR)o of Eq. (fl3|) are determined by the proton 
and neutron excess densities of Eq. (jH^ . respectively. 
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The Eq. (I12p relating the Fermi and "Coulomb" strength functions holds also within the pn- 
QRPA if the all above-mentioned conditions are fulfilled. In particular, the use of a truncated 
basis of the single-particle states within the BCS model leads to a unphysical violation of the 
isospin symmetry. In such a case the degree of violating Eq. (fT^ can be considered as a measure 
of the violation. 

The equation of motion for the GT operator Y^^^ with taking the nucleon pairing into 
consideration is somewhat modified as compared to Eqs. 0, (fTUI) . According to Eqs. (jHnjl - ljH^ 
and making use of the full basis of the single-particle states along with Eq. Q for the radial 
overlap integrals we get within the pn-QRPA (F^^"-* = '}l-Ku{XTTXv){.Q]it)^)'- 

[H,Y;f^] = {Ulr\-H + iUlr\-v- (36) 

The operator {U^^)p-h is defined by the expression (fTUI) in which both the neutron excess and 
proton densities are used according to Eq. The expression for iUi~^)p-p has the form: 



iUir\-p = %7^(A„(Pi^)+ - A,>J. (37) 



According to Eqs. (fTUjl . and (jTfjl the expression for {EWSR)i consists of two terms, the 
former, {EWSR)i ^, coinciding with Eq. (fT^ with the nucleon densities and the occupation 
numbers appropriately modified by the nucleon pairing, and the latter being due to the nucleon 
pairing only: 

{EWSRY{' = ° ^ r (38) 

Go Go 

In the SU(4)-symmetry limit one has the equality {EWSR)i = {EWSR)o, which holds also 
for double-closed-shell nuclei. 



3.3. The pn-QRPA equations. 

The system of the homogeneous equations for the forward and backward amplitudes X^^{s) = 
{s,JiJ,\ (^A^^^ |0) and Yjj^{s) = (s, J/i|A;J^|0), respectively, is usually solved to calculate the 
energies utg and the wave functions \s, Jjj) of the isobaric nucleus within the quasiboson version 
of the pn-QRPA (see, e.g., |2|). In particular, the system of the equations for the amplitudes 
follows from the equations of motion for the operators and A making use of the Hamiltonian 
(IHn jl . ljHHjl . ljn^ . Instead, we rewrite the system in equivalent terms of the elements r~'^gf{s,r) 
of the radial transition density. The elements are determined by the amplitudes X{s) and Y{s) 
as follows: 

gUs,r) = Y.it)HXAr)x.{r){s,Jl^\ {R'X). |0), z = 1,2,3,4; (39) 

where the operators P and Q are defined after Eqs. ()33|) . ()34|) . According to the defini- 
tion (jnni), the elements Qi, §2, Qs, can be called, respectively, the particle-hole, hole-particle, 
hole-hole and particle-particle components of the transition density, which can be generally 
considered as a 4-dimensional vector: {f?/}. In particular, the particle-hole strength of the 
state \s,Jfi) corresponding to a probing operator Vj^^ is determined by the element g(: 

(^j~^)s = (^TT / f?i (s, r)Vj(r) dr] (compare to Eq.(|^). The pn-QRPA system of equations 
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for the elements gf is the following (hereafter the "symmetric" approximation is only consid- 
ered) : 

1 f 

Qiis, r) = ^T.J 47*l(^' ^ = ^s)Fi!{ru r2)gi{s, ra) dndr2, (41) 

k 

F/(ri,r2) = F/(ri,r2) = 2f /^''' ~ '''\ F,'{n,r2) = FUn.r^) = -^TiGj, (42) 

rir2 

4:i(^i,^2;^) = E 0(5H)'x^(ri)x.(ri)x.(r2)x.(r2)4i(^), (43) 



A^~) -uv + ^) A^-^ - u V + 

,2 „,2 







{-} 


{+} 




'{-} 






vlul 


{-} 


{+} 




^ + 
{-} 


= Mj.fi. ( 




^^^^ 


{-} 


{+} 



{-} {+}■ 



_ _ .(-) _ , 



{-} { + }' 



^(-) „ i , - '''''' - 

^7riy,24 — "jy t^j/ 1 J- _ -r J , ^7ri.,33 ~ |_| |_|_| 5 

.(-) _V_]^_V^ (_) _.{-). .(-) _ .(-) 

where { — } = uj^""^ — E^, — E^, {+} = cu^"^ + E-„ + Ey^ u = u^^^ + /ip — /^n gives the excitation 
energy with respect to the mother-nucleus ground-state energy. All actual calculations in 
present work are performed in terms of u, whereas the experimental energy difference between 
the ground states of daughter and mother nuclei is used to represent the results in terms of 
daughter-nucleus excitation energy E^. 

According to Eq. (jUj), schematically represented as ^/ = JZk^j'iki^ — ^s)FkQk^ 4x4 
matrix (47rrJr|)~Mj^^(ri, r2; t^) is the radial part of the free two-quasiparticle propagator, 
whereas the quantities v^{s, ri) = / Fj! (ri, r2)g{{s, r2) dr2 are the elements of the radial transi- 
tion potential. To describe the excitations in /^"'"-channel one should use the same system (jlT|l 
with A^~^ superseded by A^^\ The expression for the particle-hole propagator can be 
obtained from Eqs. ()43p by the substitutions u ^ tt, n ^ p. In the case of neglecting nucleon 
pairing {Gj = 0), the system of equations (PTj) decouples and becomes determined by 

Eq. (fT^ taken in the "symmetric" approximation (it can be shown easily using the spectral 
expansion for the radial Green's functions g(^x){r,r'; e)). 

The expression for the elements of the free two-quasiparticle propagator can be obtained 
making use of the regular and anomalous single-particle Green's functions for Fermi-systems 
with nucleon pairing in analogue way as it was done in monograph to describe the Fermi- 
system response to a single-particle probing operator acting in the neutral channel. Namely, the 
matrix A can be depicted as a set of the diagrams shown in Fig.^ To illustrate this statement, 
let us consider the system of equations for the elements of the transition potential f j. The system 
follows from Eq. (^1]) and can be schematically represented as Vi = X]fc Fi^ik\^ = ^s)vk ■ 

Using expression ()43|) for the free two-quasiparticle propagator, one can obtain the sys- 
tem of the inhomogeneous equations of the pn-QRPA to calculate the strength functions ()lip 
with taking the full basis of the single-particle states into consideration for the particle-hole 
channel. The system has the form similar to that of Eqs. (|17|). (I18|) taken in the "symmetric" 
approximation: 

S?\uj) = -^ImE / Vj{r)Afl{r,r';uj)vjf{r',uj)drdr', (44) 
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V^f[r,u:) = Vj{r)5,, + ^^Y. j Ff{rr,)Afi{n,r,;u;)Vjf{r2,u;)dr,dr2. (45) 

k 

The same as in case without pairing, the residues {r^j^^)s in poles of the strength functions (see 
Eq. (1201)) coincide with the particle- hole strengths of the pn-QRPA excitations corresponding 
to the probing operators VjJ^ 



3.4. Strength functions in single-open-shell nuclei. 

We restrict the further consideration to the case of the single-open-shell nuclei. For this 
case the system of equations (ji^ for the radial effective operators vj^\r,u;) is noticeably 
simplified. Let us consider for definiteness the excitations in the /^"-channel for a nucleus in 
which only the proton pairing is realized (A„ = 0). Then, according to Eq. ()43p. the elements 
Aj~^2! ^j~i3! ^J245 ^J34 of the free two-quasiparticle propagator equal to zero and the system 
is reduced to the system of equations only for 2 radial effective operators vj^^ and vj^- 

The elements A^jJ-^^, A^^J^, A^^J^^ entering this system can be rewritten in the form which allows 
one to take the full basis of the single-particle states in the particle-hole channel for both 
subsystems. To get this form one should put E.„ = \e.j, — for E.„ > kAp {k ^ 1) and use the 
spectral expansion for the radial single-particle Green's functions in Eq. ()43|) . As a result, we 
come to the following representation for 
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^j,ii('^i,'^2;^) = wV)) nyXu{ri)Xu{r2)g{^){ri,r2]eu + u)) + (46) 



(^(i)H)^^^X^(n)x^(^2)^H(ri,r2;/i^ - - u) + 



^ {t[t){u)y ^'^Xu{ri)xAr2)xAn)xAr2) x 

..2 ^2 1 



+ 



[LJ — — Et^ + Sy LO — fip + Ej^ + £u ^ — £tv + ) 

where Yl! v-k denotes the sum, which is taken over the proton states with Et^ < kAp only. For the 
consistency of the consideration, the same truncation should be used in realization of the BCS 
model according to the second of Eqs. (p?T| . The corresponding expressions for the elements 
J 14 ^■'^cl ^j44 are! 

Aj~uirur2;uj) = -Y.' (^Mm) u^^7TxAri)X7rir2)g{u)irur2; jJ,^ - E^ - uj) + (47) 



{v)lT 



E'(*lSw)'«.AMr,)xMr.)x-,(n)x.(r.)' ' ' 



U - ^p - Et, + £y to - Hp + Et, + Eu j 

A^j'}^{ri,r2;uj) = Y' {t[n){u)) ulxn{ri)xAr2)9i,y){rur2; fi^ - E^ - uj) + (48) 



Y' {t\t)(u)y nyXu{ri)xAr2)xAri)xAr2) 



+ 



U - ftp - E^ + Ey U - flp + E^ + Ey 



The equations ()44|) - (PHj) realize a continuum version of the pn-QRPA (pn-QCRPA) on the 
truncated basis of the BCS problem in proton-open-shell nuclei. 

To go to the case when only the neutron pairing is reahzed, the equations (jiUjl - fjiHI) should 
be changed by the substitution vr z/, p ^ n, uo ^ —uj. As before, the expression for the 
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elements of the free two-quasiparticle propagator A^^^ can be obtained from the corresponding 
elements A^~^ by the substitution n ^ u, p ^ n. 



4. Calculation results 

4.1. Choice of the model parameters 

Parametrization of the isoscalar part of the mean field Uo{x) along with the values of the 
parameters used in the calculations have been described in details in j^H]- The dimensionless 
intensities fj of the Landau-Migdal forces of Eq. Q are chosen as usual: Fj = /j-300 MeVfm^. 
The value /o = /' = 1.0 of the parameter /' determining the symmetry potential according to 
Eq. ((7j) is also taken from the where the experimental nucleon separation energies have 
been satisfactorily described for closed-shell subsystems in a number of nuclei. The value of the 
dimensionless intensity /i = g' of the spin-isovector part of the Landau-Migdal forces g' = 0.8 
is chosen to reproduce within the CRPA the experimental energy of the GTR in ^°^Pb parent 
nucleus [T7] . 

The strength of the monopole particle-particle interaction Gq is chosen to reproduce the 
experimental pairing energies (actually, we identify the pairing energy with A, obtained by 
solving the BCS-model equations (jHT|) ). The summation in the second of the equations is 
limited by the interval of 7 MeV above and below the Fermi level. The same truncation is 
used in the expressions like fimjl - PHj) for the elements of the free two-quasiparticle propagator. 
Comparing the calculated value of the nucleon separation energy B'^"-'-'^ ~ — /i^ + with the 
corresponding experimental one can be considered as a test of the used version of the BCS 
model. The parameters of the BCS model along with 5^"^'^ and B^^^ are listed in the Table ^ 
Note also, that the strength of the spin-spin particle-particle interaction Gi is chosen equal 
to Go, which seems close to the realistic ones used in literature. Such a choice allows us to 
simplify the calculation control using the {EWSR)i, because the "particle-particle" part of this 
sum rule goes to zero in accordance with Eq. Another reason for the choice is the "soft" 
spin-isospin SU(4)-symmetry, which is roughly realized in nuclei (see, e.g., pSj). 



4.2. F (C) strength functions 

The F(C) strength functions have been calculated within the CRPA for the ^osp^^ parent 
nucleus according to Eqs. (jl7p . (jl8p and within the pn-QCRPA for ^^Zr and tin isotopes 
according to Eqs. (^3)) or their modifications. The "cut-off" parameter k = 3 was used in 
calculations of the radial propagators of Eqs. ()46|) - (PHj) . All above-mentioned equations have 
been taken for J = 0. The following characteristics of the F strength distribution have been 
deduced from the calculated F strength functions: the lAR energy Ejar] the mean energy 
-^/VMi? of IVMR*^~); relative Fermi strength x for 3 energy intervals: the vicinity of the 
lAR, IVMR(-), and IVMR(+); relative energy-weighted Fermi strength y for these 3 energy 
regions. As for the low-energy interval, for all nuclei in question the values of x are below 
0.1% (except for ^°Zr, where x ~0.2%). The value of {EWSR)o entering the definition of y 
has been calculated according to Eq. (jl5|) with the nucleon densities modified by the nucleon 
pairing. All the characteristics of the F strength distribution are listed in the Table |21 along 
with the experimental lAR energies. To check the selfconsistency of the model, the parameter 

^ S^~\uj)dw 



J 5^ '(a;)a;2 dui ^ 



has also been calculated. 
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4.3. GT strength functions 

The GT strength functions have been calculated within the CRPA for the ^'^^Pb parent nu- 
cleus according to Eqs. (fT7|) . ()18|) and within the pn-QCRPA for ^°Zr and tin isotopes according 
to Eqs. (j44|) ■ (j45j) or their modifications taken for J = 1. /^From the calculated GT strength 
functions the relative GT strength x and relative energy-weighted GT strength y for 4 energy 
intervals (low-energy, the vicinity of the GTR, high-energy, and the vicinity of the IVSMR(+)) 
have been deduced. The value of {EWSR)i ^ entering the definition of y has been calculated 
according to Eq. ()16p with the nucleon densities and the occupation numbers modified by the 
nucleon pairing. The value of {EWSR)i'^ defined according to Eq. (jHHj) equals zero because 
of putting Gi = Gq in our calculations. The values of the parameters x and y are listed in 
the Table El along with the calculated energies of the GTR. The values of the relative strengths 

= fs/{N — Z) of the GT doorway states calculated according to Eq. ()20|) with and with- 
out taking the nucleon pairing into consideration are shown in the Figs. |2^E] (solid and dotted 
vertical lines, respectively). 

The coupling of the GT doorway states with many-quasiparticle configurations is taken 
into consideration phenomenologically and is described in average over the energy in terms 
of the smearing parameter I{uj) simulating the mean spreading width of the doorway states. 
Following J7] we choose I {Ex) to be an universal function revealing saturation at rather high 
excitation energies: 

'(^^'> = -TTilW 

where a and B are adjustable parameters, E^ is the excitation energy in the daughter nu- 
cleus. Equation is close to the parametrization of the intensity of the optical-potential 
imaginary part obtained within the modern version of an optical model. The choice of the 
parameters a = 0.09 MeV~^ and B = 7 MeV has allowed us to describe satisfactorily the 
total widths of a number of the isovector giant resonances (including GTR) in the ^osp^ pa^-. 
ent nucleus jTHI. In this work the energy-averaged GT strength function Si~^ is calculated 
as S[~\uj) = S[~\uj + iI{Ex)/2) using the same parametrization of /, where S[~^ is the GT 
strength function calculated according to the CRPA or pn-QCRPA equations. The calculation 
results are shown in the Figs. (thin line). The calculated energy dependence in the vicin- 
ity of the GTR is approximated by the Breit-Wigner formulae to get both the GTR energy 
Egtr and width Tgtr- The values of Egtr obtained thereby and Tgtr along with the values 
calculated without taking the nucleon pairing into consideration (in brackets) are listed in the 
Table IHl in comparison with the corresponding experimental data. 

We analyze the effect of the isospin splitting according to Eqs. (j211), without taking 
the nucleon pairing into consideration because this effect turns out to be weak even for such 
a rather light parent nucleus as ^^Zi. In this case there is the only Tq component of the GT 
strength function in ^°Nb (due to the (79/2 (77/2 Ml transition in the proton subsystem). Its 
relative strength and energy are x> = 4.7 % and E^ = 12.8 MeV, respectively. It leads to a 
decrease in the relative strength and energy of the GTR equal to x> and a;>(i?> — i?<)/x< = 0.3 
MeV, respectively, as compared with the values found within the CRPA. 

5. Discussion of the results 

5.1. Summary concerned to the approach 

We have extended the CRPA approach, used previously in [T7HTU] for describing the GTR 
and the lAR in closed-shell nuclei, and have formulated a version of the pn-QCRPA approach 
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for describing the GT and F strength distributions in open-shell nuclei. The common ingredi- 
ents of both approaches are the following: the phenomenological isoscalar part of the nuclear 
mean field; the isovector part of the Landau-Migdal particle-hole interaction; the isospin self- 
consistency condition being used to calculate the symmetry potential; mean Coulomb field 
calculated in the Hartree approximation; the use of the full basis of the single-particle states 
in the particle-hole channel and, as a result, formulating the continuum versions of the RPA; 
a phenomenological description of the coupling of the particle-hole doorway states with many- 
quasiparticle configurations in terms of the mean doorway-state spreading width (the smearing 
parameter). 

The specific feature of the developed pn-QCRPA approach is the use of an isospin-invariant 
particle-particle interaction to describe both the nucleon pairing phenomenon in neutral chan- 
nels and particle-particle interaction in charge-exchange channels. The method to check the 
isospin selfconsistency, which could be violated by the use of a truncated basis of single-particle 
states in particle-particle channel within the approach, is also used. 

5.2. F (C) strength distribution 

The present results (Table ^ show at the microscopical level that the isospin is a good 
quantum number for medium-heavy mass nuclei. In particular, the JAR exhausts almost all 
F strength. The rest is mainly exhausted by the IVMR. The heavier a nucleus is, the more 
the relative F strength of the IVMR becomes. Similar results have also been obtained for the 
energy-weighted Fermi-strength distribution. 

The appropriate choice of the particle-particle interaction along with the use of the same 
truncated basis of single-particle states in the neutral and charge-exchange channels ensure the 
isospin selfconsistency of the approach (the calculated values of z are found to be less than 
0.1%). For the same reason, the effect of the nucleon pairing on the F strength distribution is 
small. 

The systematic lowing of the lAR energy compared with the experimental value is a short- 
coming of the approach, possibly caused by the absence of the full selfconsistency. 

5.3. GT strength distribution 

For all nuclei in question, the calculated GT strength function splits into three main regions 
(see Table El and Figs. I2H11)- The main part of the GT strength (70-80%) is exhausted by the 
GTR. In the case of ^°^Bi, the calculated GTR relative strength is found to be in a reasonable 
agreement with the respective experimental value Xqj^^ = (60 ± 15)% [4 . The spin-quadrupole 
part of the particle-hole interaction is taken into account only for ^°'^Bi. As a result, the cal- 
culated GTR relative strength is decreased from 68% to 66% and the low-energy GT strength 
is noticeably redistributed. In the case of ^°Nb, the difference between experimental and cal- 
culated results is somewhat larger: Xgtr = (39 ± 4)% 0, Xctr = (61 ± 8)% {x^'^p = (66l?[!)% 
up to 20 MeV in excitation energy) 0. This could be partly due to the well-known quenching 
effect which has been discussed during the last two decades. This effect is not finally established 
experimentally well enough and its discussion is out of the scope of this work. 

The calculated high-energy part of the the GT strength distribution is mainly exhausted by 
the IVSMR and its satellites. This part is relatively large for ^°'^Bi (about 19%) and decreases 
up to about 9% for ^^Nb (4.6% due to the IVSMR and 4.5% due to r> component of the GTR). 
The low-energy part contains the weakly-collectivized GT states, which can be related to those 
found in ^TJ ^1 • The nucleon pairing leads to noticeable enriching of the calculated low- 
energy part of the the GT strength distribution in open-shell nuclei. This fact is in a qualitative 
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agreement with respective experimental data (shown in Figs. |2}^- The calculated values 
xiow/xqtr 19.2% and 25.8% are found to agree reasonably with the corresponding experimental 
values 28.2% and 30.4% |7j for ^°Nb and ^°^Bi, respectively. Bearing in mind also the above- 
mentioned calculated and experimental values of xgtr, we can expect that the quenching effect 
is not so noticeable as it is usually assumed. 

The configurational splitting of the main GT state is found for some single-open-shell nuclei 
(Figs. El El). In the case of the ^°Zr parent nucleus, the splitting is caused by the direct spin-flip 
transition I/7/2 I/5/2' which becomes possible due to the proton pairing and whose energy 
is close to the GTR energy calculated without taking this transition into account. The value 
of the splitting (about 0.5 MeV) is found to be rather small. For this reason, the total GTR 
width is only slightly increased. In the case of the Sn isotopes, the strongest effect caused 
by the 1/^11/2 ~^ ^^9/2 transition is found for the ^^°Sn parent nucleus. The splitting energy 
is found to be rather large and leads to a noticeable increase of the total GTR width. The 
calculated isospin sphtting energy for ^°Nb (4.0 MeV) is found to be in good agreement with 
the respective experimental one (4.4 MeV [HI). 

In conclusion, we incorporate the BCS model in the CRPA method to formulate a version 
of the partially self-consistent pn-QCRPA approach for describing the multipole particle-hole 
strength distributions in open-shell nuclei. The approach is applied to describe the Fermi and 
Gamow- Teller strength functions in a wide excitation-energy interval for a number of single- and 
double-closed-shell nuclei. A reasonable description of available experimental data is obtained. 

Authors are grateful to Prof. M. Fujiwara and Prof. A. Faessler for thorough reading the 
manuscript and many valuable remarks. V.A.R. would like to thank the Graduiertenkolleg 
"Hadronen in Vakuum, in Kernen und Sternen" (GRK683) for supporting his stay in Tubingen 
and Prof. A. Faessler for hospitality. 
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Table 1: Calculated values of the pairing gap A/? along with the calculated and experimental 
proton and neutron separation energies -B™^'^' The monopole pairing strength Gq is also 
given 



Nucleus 


Go- A 

MeV 


A/3, 
MeV 


ID calc 

MeV 


MeV 


MeV 


MeV 


QHr? 

^"Zr 


14.0 


1.15 


11.5 


8.2 


11.98 


8.36 


ii^Sn 


11.5 


1.3 


10.9 


7.5 


10.79 


7.56 


ii^Sn 


11.5 


1.4 


10.7 


8.3 


10.30 


8.48 


ii^Sn 


11.5 


1.3 


9.9 


9.1 


9.56 


9.28 


ii^Sn 


11.5 


1.4 


9.7 


9.9 


9.33 


10.00 




11.5 


1.4 


9.4 


10.7 


9.11 


10.69 


i22Sn 


11.5 


1.4 


9.1 


11.5 


8.81 


11.4 


i24Sn 


11.5 


1.2 


8.7 


12.3 


8.49 


12.10 



Table 2: Calculated characteristics of the Fermi strength distribution along with respective 
experimental data IH El (see the text for details) 



Parent 
nucleus 


X, % 


y, % 


EjAR , MeV 


ElVMR , 

MeV 


lAR 


high 


(+) 


lAR 


high 


(+) 


calc. 


exp. 


902r 


98.6 


1.5 


0.5 


94.2 


4.2 


0.9 


4.6 


5.0 


27.1 


i^^Sn 


97.8 


2.4 


0.6 


92.9 


5.6 


0.7 


5.6 


6.2 


25.0 


ii^Sn 


98.0 


2.4 


0.5 


93.5 


5.5 


0.7 


6.8 


7.3 


25.3 




97.3 


2.4 


0.4 


93.2 


5.4 


0.5 


7.8 


8.4 


25.5 


ii^Sn 


97.7 


2.4 


0.3 


93.9 


5.3 


0.5 


8.8 


9.3 


26.0 




97.6 


2.4 


0.3 


93.9 


5.2 


0.4 


9.7 


10.2 


26.4 


i22Sn 


97.5 


2.4 


0.2 


94.0 


5.2 


0.3 


10.7 


11.2 


26.8 




96.8 


2.5 


0.2 


93.4 


5.2 


0.3 


11.7 


12.2 


27.3 


208pjg 


94.4 


3.1 


0.2 


89.8 


6.2 


0.2 


14.4 


15.1 


32.5 
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Table 3: Calculated characteristics of the Gamow- Teller strength distribution along with re- 
spective experimental data O IH |7j (see the text for details) . 



Parent 
nucleus 


X, % 


y, % 


Egtr, MeV 


Vgtr. MeV 


low 


Lr i K 


high 


(+) 


low 


(or i K 


high 


(+) 


calc. 


exp. 


calc. 


exp. 




1 /I s 
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Figure 1: Diagrams corresponding to elements A^^J^l^ of the free two-quasiparticle propogator 
of Eq. PHj) . Upper (lower) lines correspond to propagation of proton (neutron) quasiparticles. 
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Figure 2: The relative GT strengths x in ^°Nb, ^^^~^^^Sb calculated within the pn-QCRPA (solid 
vertical line) and CRPA (dotted vertical lines) in comparison with the respective experimental 
data 13 El (dashed vertical lines). The smeared pn-QCRPA GT-strength distribution S is also 
shown (thin solid lines) 
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Figure 3: The same as in the Fig. |21 but for the ^^^Sb (the experimental data are taken 
from [^) 
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Figure 4: The relative GT strengths x in ^''^Bi calculated within the CRPA (solid vertical lines) 
along with the respective smeared distribution S (thin solid line). The experimental data are 
taken from [7j (dashed vertical lines) 
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